International Journal of Engineering and Technical Research (IJETR) 
ISSN: 2321-0869 (O) 2454-4698 (P), Volume-5, Issue-2, June 2016 


Influence of radiation on MHD free convective flow 
of a Williamson fluid in a vertical channel 
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Abstract — In this paper, we studied the effects of MHD and 
radiation on fully developed free convection flow of a 
Williamson fluid in a vertical channel. The governing non-linear 
equations are solved for the velocity field and temperature field 
using the perturbation technique. The effects of various 
emerging parameters on the velocity field, temperature field, 
skin friction and Nusselt number are studied through graphs 
and tables in detail. 

Index Terms — Williamson fluid, MHD, Nusselt number, 
radiative heat flux, Hartmann number 


I. INTRODUCTION 

The insufficiency of the classical Navier-Stokes theory to 
describe Theologically complex fluids such as blood, paints 
and polymeric solutions, has led to the development of several 
theories of non-Newtonian fluids. Non-Newtonian fluid 
flows are encountered in a wide range of engineering 
applications. Hot rolling, extrusion of plastics, flow in journal 
bearings, lubrication, and flow in a shock absorber are some 
typical examples. The increase of these applications in the 
past few decades has urged scientists and engineers to provide 
mathematical models for non-Newtonian fluids. The 
nonlinearity between stress and deformation rate for this kind 
of fluids makes it, in general, impossible to obtain a simple 
mathematical model as in the case for Newtonian fluids. This 
difficulty has lead researchers to investigate relatively simple 
non-Newtonian fluid models. A detailed discussion on second 
and third-order fluids can be found in the study of Dunn and 
Rajagopal (1995). Hady and Gorla (1998) have investigated 
the effect of uniform suction or injection on flow and heat 
transfer from a continuous surface in a parallel free stream of 
viscoelastic second order fluid. The flow of third grade fluid 
between heated parallel plates was studied by Akyildiz 
(2001). Chamka et al. (2002) have discussed the fully 
developed free connective flow of micropolar fluid between 
two vertical parallel plates analytically. Siddiqui et al. (2010) 
have analyzed the flow of a third grade non-Newtonian fluid 
between two parallel plates separated by a finite gap by using 
the Adomian decomposition method. Williamson fluid is 
characterized as a non-Newtonian fluid with shear thinning 
property, i.e., viscosity decreases with increasing rate of shear 
stress (Dapra and Scarpi, 2007). 

The study of free convective flow in a channel of an 
electrically conducting viscous fluid under the action of a 
transversely applied magnetic field has immediate 
applications in many devices such as magnetohydrodynamic 
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(MHD) power generators, MHD pumps, accelerators, 
aerodynamics heating, electrostatic precipitation, polymer 
technology, petroleum industry, purification of crude oil and 
fluid droplets sprays. Attia and Kotb (1996) have studied the 
MHD flow between two parallel plates with heat transfer. 
Free convective flow of micropolar fluid between two parallel 
porous vertical plates with the effect of magnetic field was 
investigated by Bhargava et al. (2003). Hayat et al. (2004) 
have analyzed the Hall effects on the unsteady hydromagnetic 
oscillatory flow of a second grade fluid. The unsteady flow of 
a dusty conducting fluid between parallel porous plates was 
studied by Hazeem attia (2005). Sanyal and Adhikari (2006) 
have investigated the effects of radiation on MHD fluid flow 
in vertical channel. The transient MHD couette flow of a 
Casson fluid between parallel plates with heat transfer was 
discussed by Attia and Ahmed (2010). Arunakumari et al. 
(2012) have investigated the fully developed free convection 
flow of a Williamson fluid in a vertical channel under the 
effect of magnetic field. 

In view of these, we studied the effects of MHD and 
radiation on fully developed free convection flow of a 
Williamson fluid in a vertical channel. The governing 
non-linear equations are solved for the velocity field and 
temperature field using the perturbation technique. The 
effects of various emerging parameters on the velocity field, 
temperature field, skin friction and Nusselt number are 
studied through graphs and tables in detail. 


II. MATHEMATICAL FORMULATION 




V 


Fig. 1 The physical model 
We consider the laminar free convection flow of a 
Williamson fluid between two plates at distance h a apart 
with radiative heat transfer, as shown in Fig. 1.We choose 
co-ordinates system, with X - axis parallel to the flow while 
Y - axis is normal to the flow. A uniform magnetic field B {) 
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is applied perpendicular to the fluid flow. The flow is 
assumed steady and fully developed, i.e., the transverse 
velocity is zero. It is also assumed that the walls are heated 
uniformly but their temperatures may be different resulting in 
asymmetric heating situation under these assumptions the 
equations that describe the physical situation are 


d 2 u ,, d 

%-7T + % r -r 
dy ay 


\2 


- crBlu + pgj3(T-T 0 ) = 0 


( 1 ) 


k <?T 1 dq _ 

dy 2 c dy 


> ~p (2) 

where p is the density, g is the acceleration due to gravity, 

P coefficient of thermal expansion, T is the viscoelastic 

parameter, C is the electrical conductivity , C p is the 

specific heat at constant pressure, k is the thermal 
conductivity and q is the radiative heat flux. Following 
Cogley et al. (1968), it is assumed that the fluid is optically 
thin with a relatively low density and the radiative heat flux is 
given by 

^L=4a 2 (T 0 -T) (3) 

Equations here CL is the mean radiation absorption 
coefficient. 

Subject to the boundary conditions 
w(0) = 0, r(0) = 2;, u(h) = 0, T(h) = T 2 (4) 
Introducing the following non-dimensional variables 
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into Eqs. (1) and (2), we get (after dropping the bars) 
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where M — B h I _ is the Hartmann number, 
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Gr — ------— is the Grashof number, We =- 


is the Wesseinberg number, N = 


2 ah 

pk 


is the radiation 


TJh 

parameter and R e = ___ is the Reynolds number, 
v 

The corresponding dimensionless boundary conditions 

w(0) = 0, 0(0) = r T , w(l) = 0, 6>(1) = 1 

( 8 ) 


III. SOLUTION 

Eq. (6) is non-linear and it is difficult to get a closed form 
solution. However for vanishing Wesseinberg number We , 
the boundary value problem is agreeable to an easy analytical 
solution. In this case the equation becomes linear and can be 
solved. Nevertheless, small T suggests the use of 


perturbation technique to solve the non-linear problem. 
Accordingly, we write 

u = u 0 + Weu x (9) 

and 0 = 0^+ WeO x (10) 

Substituting equations (9) and (10) into Eqs. (6) and (7) and 
boundary conditions (8) and then equating the like powers of 
We , we obtain 

A. Zeroth-order system (We°) 
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i (0) = u o (l) = 0 ’ #0 (0) = r T , # 0 (l) = 1 (13) 

B. First-order system (We) 
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Together with boundary conditions 

m 1 (0) = m 1 (1) = 0. 6>(0) = 0, 6>(1) = 0 


(16) 


C. Zeroth-order solution 

Solving Eqs. (3.3) and (3.4) using the boundary conditions 
(3.8), we get 

0 O = A l cos Ny + A 1 sin Ny (17) 

rz r 
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D. First-order solution 

Solving Eq. (3.7) subject to the boundary conditions in Eq. 
(3.8), we get 

0 X = 0 (19) 

Substituting the Eqs. (3.10) and (3.11) into the Eq. (3.6) and 
then solving the resulting equation with the corresponding 
conditions, we get 
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Finally, the perturbation solutions up to first order for 0 
and U are given by 

6 = 0 0 + We6 x = A cos Vy + A sin Ay ( 21 ) 
and U = Uq 4- Weu x 


( 22 ) 

The skin friction at the plate y = 0 of the channel is given 
by 


T 


du 



diir, du, 
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The rate of heat transfer coefficient in terms of Nusselt 
number Nu at the plate y = 0 of the channel is given by 


A, 50 

Nu = — 
dy 


= N^ 


y=0 


(24) 


IV. RESULTS AND DISCUSSIONS 


In order to perceive the effects of various physical 
parameters like Gr , Re , We , M , N and r T on the 

velocity profile u we have plotted the Figs. 2-7. From Fig. 2, 
it is observed that the velocity U first decreases and then 
increases with increasing Wesseinberg number We. From 
Fig. 3, it is seen that the velocity U decreases with increasing 
Hartmann number M . From Fig. 4, it is found that the 
velocity U increases with an increase in radiation parameter 
N . From Fig. 5, it is observed that the velocity U increases 
with increasing Grashof number Gr . From Fig. 6, it is noted 
that, the velocity U decreases with an increase in Reynolds 
number Re. From Fig. 7, it is found that, the velocity U 
increases with increasing V T . 

Table-1 depicts the effects of various physical parameters 
like Gr , Re , We, M , N and r T on the skin friction T 
at the plate y — 0 . It is found that, the skin friction T 
increases with increasing radiation parameter N , Grashof 
number Gr and r T , while it decreases with increasing 

Wesseinberg number We , Hartmann number M and 
Reynolds number Re. 

Fig. 8 shows the effect of wall temperature parameter V T 
on 0. It is observed that, the temperature 6 increases with 
an increase in V T . 

Fig. 9 depicts the effect of radiation parameter N on 6. it 
is observed that, the temperature 6 increases with an 
increase in V T . 

Table-2 shows the effect of V T on the Nusselt number 
Nu at the plate y — 0. It is found that, the Nusselt number 
increases with increasing radiation parameter N, while it 
decreases with increasing r T . 
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Fig. 2. Effect of Weisseiiberg number We on // for Gr = 1, N = 1, M = 1, 
r T = 0.5 and Re = 1. 



Fig. 3. Effect of Hartmann number M on u for Gr = 1, N = 1. We = 0.1, 
r T =0.5 and Re = l. 



Fig. 4. Effect of radiation parameter N on it for Gr = 1. We = 0.1. M = 1, 
7 j =0.5 and Re = 1. 



Fig. 5. Effect of Grashof number Gr on it for N = 1. We = 0.1. M — 1, 
r T = 0.5 and Re = 1. 



Fig. 6. Effect of Reynolds number Re on it for Gr = 1. We = 0.1. M = 1, 
7 j =0.5 and N = 1. 



Fig. 7. Effect of r T on n for Gr = 1, We = 0.1. M = 1, tV = 1 and Re = 1. 



Fig. 8. Effect of wall temperature parameter r T on 0 for N = \ . 



y 

Fig. 9. Effect of radiation para meter N on d for r T =0.5 . 
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Til Ible — 1: The Skill friction r at the plate y = 0 
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Table — 2: The Nusselt number Nil at the plate y = 0 


N 

V 7 

Nit 

1 

0.5 

0.8673 

2 

0.5 

2.6572 

1 

0.8 

0.6747 


V. CONCLUSION 

In this chapter, we studied the MHD fully developed free 
convection flow of a Williamson fluid in a vertical channel 
with radiative heat transfer. The governing non-linear 
equations are solved for the velocity field and temperature 
field using the perturbation technique. It is found that, the 

velocity increases with increasing N,Gr and V T , while it 
decreases with increasing We , M and Re ; the skin friction 
T increases with increasing a , Gr and r T , while it 
decreases with increasing We , M and Re ; the temperature 
0 increases with an increase in V T and the Nusselt number 
decreases with increasing r T . 
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